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ABSTRACT 

High-resolution submillimetre (submm) imaging of the high-redshift radio galaxy (HzRG), 
4C 60.07, at z = 3.8, has revealed two dusty components of roughly equal integrated flux. 
Spitzer imaging shows that one of these components ('B') is coincident with an extremely 
red active galactic nucleus (AGN), offset by ~4arcsec (~30kpc) from the HzRG core. The 
other submm component ('A') - resolved by our synthesised beam and devoid of emission 
at 3.6-8.0-/im - lies between 'B' and the HzRG core. Since the radio galaxy was discovered 
via its extremely young, steep-spectrum radio lobes and the creation of these lobes was likely 
triggered by the interaction, we argue that we are witnessing an early-stage merger, prior to 
its eventual equilibrium state. The interaction is between the host galaxy of an actively-fueled 
black hole (BH), and a gas-rich starburst/AGN ('B') marked by the compact submm compo- 
nent and coincident with broad CO(4-3) emission. The second submm component ('A') is a 
plume of cold, dusty gas, associated with a narrow (^150km s _1 ) CO feature, and may repre- 
sent a short-lived tidal structure. It has been claimed that HzRGs and submm-selected galaxies 
(SMGs) differ only in the activity of their AGNs, but such complex submm morphologies are 
seen only rarely amongst SMGs, which are usually older, more relaxed systems. Our study 
has important implications: where a galaxy's gas reservoir is not aligned with its central BH, 
CO may be an unreliable probe of dynamical mass, affecting work on the co-assembly of BHs 
and host spheroids. Our data support the picture wherein close binary AGN are induced by 
mergers. They also raise the possibility that some supposedly jet-induced starbursts may have 
formed co-evally (yet independently of) the radio jets, both triggered by the same interaction. 
Finally, we note that the HzRG host would have gone unnoticed without its jets and its com- 
panion, so there may be many other unseen BHs at high redshift, lost in the sea of ~ 5 x 10 8 
similarly bright IRAC sources - sufficiently massive to drive a >10 27 -W radio source, yet 
practically invisible unless actively fueled. 

Key words: galaxies: starburst - galaxies: formation - techniques: interferometric - instru- 
mentation: interferometers - submillimetre 



1 INTRODUCTION 




Radio galaxies are believed to host actively-fueled, spinning 
BHs which power their immense radio luminosities and fash- 



extreme youth (^lOMyr - iBlundell & Rawlingj [l999l) . Nowa- 




1 USS sources, with a < — 1 where S„ oc i/ a , were often found to la ck 
identifications on photographic plates (e.g. Hel ens, Milev & Willis! 1979T) . 



examples of the phenomenon (HzRGs) have often been iden- 
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days, a wealth of unequivocal evidence also links HzRGs with 
very massive galaxies in the early Universe. Near- and mid- 
infrared (-IR) observations show that HzRGs are associated 
with the most massive stellar populations at any given redshift 
( iBest. Longair & Rottgerind 1 19981 : ISevmour et"atl |2007|) . Direct 
evidence of vast reservoirs of atomic and molecular hydrogen 
has also been established, via observations of strong Hi absorp- 
tion against luminous, morphologically complex Lya halos - of- 
ten extending 100-200 kpc fro m the central radio galaxy — and 
via detections of CO and dust (Hippelein & Meisenheimer 1993; 



Papadopoulos et all 200d: lArchibald et alj 


2001; Reuland et alj 


2001 2004: De Breuck et al.l 12001 2005: 


iKlamer et al.l 12005: 


Villar-Martm et al. 2006). 



If HzRGs pinpoint the most massive galaxies out to the highest 
redshifts then we expect these systems to signpost the most mas- 
sive peaks in the primordial Universe and to evolve into today's 
brightest galaxies - cD-type ellipticalfl In accord with these ex- 
pectations, the environments of HzRGs are over-dense in a rich 
variety of galaxy types, including SMGs. The first e xample, still 
the most striking, is the field of 4C41.17 at z = 3 .8 dlvison et al.1 
l200Ct iGreve et al.ll2007l) , and lStevens et alj d2003h presented low- 
resolution submm imaging of a further six HzRGs. Some of these 
displayed evidence of extended (/^ 5-20 arcsec, ^ 30-150 kpc) dust 
emission (and UV emission - lHatch et alj|2008l) . suggesting that 
starbursts in HzRGs differ from the compact events seen in lo- 
cal ultraluminous IR galaxies (ULIRGs) and, indeed, in the gen- 
eral 2 ~ 1 — 3 SMG population as probed by high-resolution 
C O and submm-radio c ontinuum imaging ( £0.5 arcsec, or £4 kpc 
IChapman et al.ll2004l; iTacconi et alj |2006, 2008; Big gs & Ivisor] 



120081; I Younger et al.l20 08b), i.e. that the mechanism for the forma- 



tion of these very massive galaxies may be fu ndamentally diffe rent 
to that of >L* galaxies forming from SMGs dSmail et alj2004h . 

Probing the size and morphology of the rest-frame far-IR 
emission in HzRGs is thus a key piece in the puzzle of galaxy 
formation. With the ~14-arcsec (FWHM) spatial resolution usu- 
ally available, we have been unable to distinguish between entirely 
different modes for the formation of the most massive galaxies. 
One could imagine HzRG host galaxies assembling via a mas- 
sive, widespread burst with a low overall star-formation density, 
or within multiple gas-rich sub-components destined to merge. Ev- 
idence from mm interferometric imaging has been enlightening, 
though the need for imaging on a range of spatial scales is often 
evident: iDe Breuck et al . (2005) found evidence for two gas com- 
ponents in the 4C41.17 system whilst IPapadopoulos et al . (2000) 
argued for a galaxy-wide starburst in another HzRG, 4C 60.07, with 
gas being consumed on a scale of ~30kpc, and for a gas reservoir 
more commensurate with local ULIRGs in another, 8C 1909+722. 

4C 60.07 is a powerful double-lobed FR II radio galaxy - one 
of t he brightest of the H zRGs observed by lArchibald et alj j200lT) 
andlReuland et aQ d2004l) in the submm waveband. The radio spec- 
trum of 4C 60.07 is unremarkable at low frequencies, with a spec- 
tral index, a = —0.9 between 38 and 178 MHz (where S v oc u a ), 
but its spectrum steepens at higher frequencies, reaching a = —1.6 
between 4.9 and 15 GHz, hence its selection as an USS radio emit- 
ter and its eventual identi fication with an J?y ega ~ 23.2 galaxy 
at 2 = 3.7887 ± 0007 dChambers et"ai]| 19961: iRottgering et"atl 
ll997l : lReuland et alj2007h . 



2 Defined as being centrally located in clusters and very much larger and 
brighter than other galaxies in the cluster, the c being analo gous to the no- 
tation for supergiant stars in stellar spectroscopy (Matthews et al. 1964). 



In this paper we present submm interfe rometry of 4C 60.07, 
utilising the Submillimetre Arra>0 (SMA - IHo et all |2004|) at a 
wavelength of 890 /im with a spatial resolution of ~2 arcsec and 
high astrometric precision - the first high-resolution submm con- 
tinuum image of an unmistakably massive galaxy at 2 ~ 4, sensi- 
tive to dust on the scales predicted by existing single-dish imaging. 
The primary goal was to determine the size and nature of the emis- 
sion. Is it diffuse, or is it due to multiple discrete, compact objects 
whose total flux density can account for the high Ssso^m seen by 
SCUBA? 

In the next section we describe an extensive set of observations 
and the associated analysis. Presentation of the reduced images fol- 
lows in Sj3] with our interpretation of those data and discussion of 
their implications in Sj4]and !j5] 



2 OBSERVATIONS 
2.1 SMA 

Data were obtained on Mauna Kea, Hawaii, during 2007 March 
and October using the SMA's "compact" configuration (baselines 
as long as 113 m during March and 70 m during October - see 
Fig. [l) and the 345-GHz receivers, with 2 GHz of bandwidth 
in each of the upper and lower sidebands, as outlined in Ta- 
ble0 The pointing centre was R.A. 05:12:54.80, Dec. +60:30:51.7 
J2000 and observations of 4C 60.07 were interspersed with ~5- 
min scans of the phase calibrator, 4C50. 11, a 2.3- Jy source lo- 
cated 14° from the target. The October data also included obser- 
vations of 8C 0716+714 (1 Jy, 17°) - to calibrate those data taken 
when the target was setting - and NVSS J044923+633208 (0.3 Jy, 
4°) to check the quality of the phase referencing. The position of 
NVSS J044923+633208 was found to be accurate to £0.1 arcsec. 

All calibration was carried out using MIR, a set of IDL-based 
routines specifically developed for data from the SMA. The main 
steps consisted of flagging noisy or corrupted data, converting the 
data to units of flux density using antenna system temperatures, de- 
termining the bandpass (with observations of 3C 454.3 or 3C 279) 
and, finally, calibrating the antenna complex gains (phase and am- 
plitude) and interpolating the calibrator solutions onto the target. 
The absolute flux scale is estimated to be accurate to 10 per cent. 
The data were then read into the Astronomical Image Processing 
Software (AZVS) package where they were imaged. 

The duration of our October track was 7.35 hr, with around 
1 mm of precipitable water vapour (pwv). The resulting noise level 
was 1.4 mjy beam" 1 , with a [2.4 x 2.2] arcsec FWHM synthesised 
beam for naturally weighted uv data. The SMA Beam Calcula- 
tor and Sensitivity Estimatoi0 had predicted a synthesised beam 
profile of [2.2 x 2.1] arcsec FWHM and an r.m.s. noise level of 
1.41 mJy beam - , for a 7.5-hr track with a 20° minimum elevation 
and 1 mm of pwv, and so appears to function adequately. 

We experimented with several permutations of data. All com- 
binations resulted in images with the same basic characteris- 
tics. Although combining all the data yielded the lowest noise, 



3 The Submillimeter Array is a joint project between the Smithsonian As- 
trophysical Observatory and the Academia Sinica Institute of Astronomy 
and Astrophysics and is funded by the Smithsonian Institution and the 
Academia Sinica. 

4 Additional data were taken in poor weather conditions during 2007 
February 20 - see TablefT]- but were not used here. 

|http : / /smal . sma. hawaii . edu/beamcalc . html| 
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Table 1. Observations of 4C 60.07 with the SMA. 



Date 


Track length 


Conditions 


Comment 


(2007) 


/min 


T225GHz 




February 20 


216 


0.167 ±0.012 


Unsuitable 


March 02 


257 


0.063 ± 0.012 


Adequate 


March 03 


268 


0.108 ± 0.008 


Marginal 


October 23 


441 


0.053 ± 0.004 


Excellent 



1.19mJybeam -1 , we opted to use only those data taken in gen- 
uine 'submm weather' (T225GHZ < 0.1), on October 23 and 
March 02, since experience suggests only these produce adequate 
representations of the submm sky. The resulting noise level was 
1.38 mjy beam - 1 , before CLEANing, with a synthesised beam mea- 
suring [2.25 x 2.07] arcsec FWHM with the major axis at position 
angle (PA), 78°. Employing 100 interations of the CLEAN algo- 
rithm, with a gain of 0. 1, in a region of radius 45 arcsec centred on 
the radio galaxy (nine times the area contained within the SMA's 
primary beam), resulted in a 20 per cent reduction in the measured 
noise. The source structure and flux density were unaffected. 



2.2 Spitzer 

4C 60.07 was observed with Spitze^ dWerner et ai]|2004l) in all the 
broadband filters a vailable in the Mu ltiband Imaging Photometer 
■ Rieke et al]|2004l) and Infrared Array Camera 



for Spitzer (MIPS 



(IRAC - iFazio et al.ll2004l) bands, 3.6-8.0 and 24-160 pm, respec- 
tively, in 2004 October, i.e. durin g the cold portion of the mission 
(programme identification: 3329 - ISevmour et alj2007b . 

The IRAC data consist of a short sequence of four dithered 30- 
s frames centred on 4C 60.07 in each of the four bands. We elimi- 
nated scattered light, column pulldown, residual images and multi- 
plexer bleed from the individual exposures by hand, then mosaiced 
these cosmetically superior frames using IRACproc dSchuster et al.l 
2006) with 0.86-arcsec pixels (half the area of the native IRAC pix- 
els). 

The Infrared Spectrograph (IRS - iHouck et al]2004l) also took 
two 61-s exposures at 16 /im in the 4C 60.07 field, using 'peak-up 
imaging' mode, prior to it being offe r ed to the community. We use 
the image reduced bv lSevmour et al.l d2007b here. 

The reduction of the 24- data was based on version S 16. 1 .0 
of the MIPS pipeline. We used an object-masked median stack of 
all the MIPS exposures of 4C 60.07 to compensate for structure 
arising from residual images. The remaining structure in the back- 
grounds was eliminated using IMSURFIT in IRAF and the resulting 
mosaics were created with 1.2-arcsec pixels. 

The IRAC, IRS and MIPS data are aligned to less than 
0.5 arcsec, conservatively. 

For IRAC and MIPS, photometric measurements were made 
using 4.9- and 7.0-arcsec-diameter circular apertures, respectively, 
correcting to total flux via the standard aperture correctiorQ. The 
flux did not vary significantly for apertures ranging from 6 to 








& 




. 



Figure 1. Top: Coverage in the uv plane of SMA data acquired on 2007 
March 02 and October 23. Bottom: Greyscale representation of the dirty 
beam, with contours at —10, —5, 5, 10, 15, 20, 40, 60 and 80 per cent of 
the peak. 



9 arcsec. At all wavelengths, apertures were centred on the 4.5-^im 
centroid. 



2.3 Optical and near-IR observations 

An /-band image was obtained during 2000 November 01 using 
the 4.2-m William Herschel Telescope 's (WHlJ3 Mosaic Cam- 
era, which co mprises two E EV 4096 x 2048 detectors with 0.236- 
arcsec pixels (Tulloch 2000). The 3-cr limiting magnitude in a 4.9- 

arcsec-diameter aperture is Ivega = 24.8; 

We also utilise the K' imaging of van Breugel et al. < 1998b . 
a 1-hr exposure taken using NIRC dMatthews & Soiferll 1994b on 
the 10-m Keck I telescope during 0.65-arcsec seeing, with a 3-a 
limiting magnitude in a 4.9-arcsec-diameter aperture of Ky cga = 
22.0. 



6 This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. Support for this work 
was provided by NASA through an award issued by JPL/Caltech. 

7 IRAC Data Handbook; http://ssc.spitzer.caltech.edu/irac/dh/ 



8 Based on observations made with the WHT operated on the island of La 
Palma by the Isaac Newton Group in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofisica de Canarias 
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2.4 Radio observations 

Pseudo-continuum data were obtained at 1.4 GHz during 2001 Jan- 
uary using the NRAO's Very Large Array (VLA) in its A con- 
figuration (programmes AD432 and AI84). The raw data were 
flagge d, calibrated and im aged following the technique described 
by Big gs & Ivisonl J2006h and the resulting images have a synthe- 
sised beam size of [1.4 x 1.1] arcsec at a PA of 26°. 

We obtained still higher-resolution continuum images, with 
100-MHz bandwidths centred at 4.7 and 8.2 GHz (programme 
AC374) and synthesised beam sizes of ~0.5 and ~0.3 arcsec 
FWHM from the National Radio Astronomy Observatory'^] 
(NRAO's) Data Archive SystenQ 

We al s o red uced the 24-GHz continuum data taken by 
iGreve et"al] d2004l) to explore CO J = 1 — emission from 
4C 60.07 for the velocity range of t he broad CO(4-3) component 
found by IPapadopoulos et al. (200fj): +455 ± 312 km s" 1 . These 
were obtained between 2001 October and 2003 March and include 
five, two and two tracks of data in the D, C and DnC configu- 
rations, respectively (programmes AI88, AI93 and AI104). Com- 
bining these data with a natural weighting scheme (ROBUST = 5) 
yielded a [2.6 x 2.4]-arcsec synthesised beam (FWHM), at PA 102°, 
and a = 11.9 /^Jy beam" 1 after 100 iterations of the CLEAN algo- 
rithm in AXVS IMAGR with a gain of 0.1 - mo re than twice as 
deep as the maps presented bv lGreve et alj J2004l) who used only 
the highest resolution data. 

We used several methods to quan tify the CO J = 1 — emis- 
sion. First, we aped the approach of IGreve et al. I d2004h : we im- 
aged the two 50-MHz-wide intermediate frequencies (IFs), repre- 
senting 'continuum only' (IF1, a = 15.7 (iJy beam -1 ) and 'con- 
tinuum plus line' (IF2, a — 15.9 /xJy beam" 1 ), then subtracted 
IF1 CLEAN components from the uv database for IF2, after check- 
ing the position and integrated flux density of the brightest lobe in 
IF1 and IF2. The dirty image made using the continuum-subtracted 
data has a — 16.0 /iJy beam" 1 . For our second method, we sub- 
tracted a dirty image of IF1 from one of IF2, yielding a ~ \/2x 
noisier map than the first method (a — 23.5 ^Jy beam ). For 
our third method, we repeated the second method for each of the 
nine observational epochs, tapering the data to produce a ~5-arcsec 
(FWHM) beam. The noise-weighted combination of these nine im- 
ages yielded a = 30.0 /^Jy beam" 1 . The three methods gaveacon- 
sistent picture in which the peak flux density due to broad-line CO 
J = 1 — emission centred on any of the prominent IR or ra- 
dio components can be no greater than 0.05 mjy beam" 1 , and no 
greater than 0. 12 mJy in total (both 3 a). 



2.5 Submm observations 

Observatio ns using the Submm Common-User Bolometer Array 
(SCUBA - iHolland et aill 19991) on the 15-m James Clerk Maxwell 
Telescope (JCMT) at a wavelength of 850 /im were described by 
IStevens et alj fc003h . We have returned to those data, and the 450- 
/j,m data obtained si multaneously, reducing them afresh with the 
software described bv llvisor] fcOOrJ) , though without the shift-and- 
add technique since 4C 60.07 was not detected at high signal-to- 
noise (S/N) in individual scans. The resulting images, smoothed 
with 7-arcsec Gaussians to yield point spread functions with FWHM 



9 The National Radio Astronomy Observatory is operated by Associated 
Universities Inc., under a cooperative agreement with the National Science 
Foundation. 

lf) https://archive.nrao.edu 



Table 2. Submm and radio positions in the 4C 60.07 system. 



• • 

Charac tens tic 


Value 


Comment 


R.A. (J2000) 


05:12:55.13 ±0.04 


NE SMA component 'A' 


Dec. (J2000) 


+60:30:49.8 ±0.3 




Peak" S 8 90/im 


6.34 ± 1.36 mjy 




Total" 5 890Mm 


8.72 ± 2.93 mJy 




R.A. (J2000) 


05:12:54.71 ±0.03 


SW SMA component 'B' 


Dec. (J2000) 


+60:30:48.8 ±0.2 




Peak" S890/jm 


7.86 ± 1.39 mjy 




Total" S890Aim 


7.83 ± 2.40 mjy 




R.A. (J2000) 


05:12:55.147 ±0.003 


Radio core, 4.7 GHz 


Dec. (J2000) 


+60:30:51.01 ±0.02 




R.A. (J2000) 


05:12:55.154 ±0.004 


Radio core, 8.2 GHz 


Dec. (J2000) 


+60:30:51.03 ±0.03 





"Measured using a twin-Gaussian fit to the dirty image. 

of ~10 and ~15 arcsec at 450 and 850 /im, are shown in Fig. 
No positional corrections have been applied. Another faint SMC 
lies to the north, with its reference beams visible, well outside the 
primary beam of the SMA's 6-m antennas at 345 GHz (30 arcsec, 
FWHM). 



3 RESULTS 

3.1 Submm morphology and flux density 

As stated earlier, 4C 60.07 is extremely bright at submm wave- 
lengths, with Sgsu/jm = 23.8 ± 3.5 mjy, as measured via fully- 
sampled jiggle maps with SCUBA (Fig. [2}- Using SCUBA's pho- 
tometry mode, which was insensitive to emission larger than 
the telescope beam, the measured flux density was 17.11 ± 
1.33 mjy beam" 1 - th e bri ghtest of 47 radio gal axies observed by 
lArchibald etHl feOOlh and lReuland et alj J2004h . 

Our new, high-resolution submm interferometry, shown 
in Fig. [3] reveals unambiguously that the emission first de- 
tected using SCUBA is spread over a region larger than the 
<500pc (<0.1 arcsec at z = 3.8) usua lly encountered in lo- 
cal ULIRGs dPownes & Solomon! 1 19981) . or eve n the £4kpc 
(<0.5 arcsec) usually encountered in SMGs (e.g. lYounger et alj 
l2007l l2008al ; iBiggs & Iyisor] 120081 cf SMM J094 303+4700 and 
SMM J123707+6214 in iTacconi et alj [200(1 120081) . Two submm 
components are revealed, hereafter 'A' and 'B'. Component 'A' is 
the easternmost and northernmost (see Table [2}, lying 1.2 arcsec 
south of the flat-spectrum radio core that has hitherto been pre- 
sumed to host the AGN driving 4C 60.07's prodigious radio lumi- 
nosity. A' and 'B' lie 3.3 arcsec apart (24kpc in the plane of the 
sky) and have similar integrated flux densities (S^o^m = 8.7 and 
7.8 mjy). Hereafter, A', 'B' and the radio core will be the principal 
components discussed. 

Component A' is marginally resolved by our ~2-arcsec syn- 
thesised beam, with a size reported by JMFIT in AXVS of [2.2 x 



0.4] arcsec at PA = 78°; its peak S^o^m is lower than that of 'B', 
but its integrated S^oo^m is marginally higher. The total Sago^m °f 



Labelled '1' in IStevens et ai] 120031) . which contains spurious coordi- 
nates. The object should be referred to as SMMJ051254.2+603119. It has 
no counterpart at 24 fim and may lie at a similarly high redshift to 4C 60.07. 
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Figure 2. Greyscale representation of 4C 60.07 observed using SCUBA 
(i |2.5> at 450 fj,m, superimposed with contours of the 8 50- fim im- 
age (shown in its original form by I Stevens et al ] I2003I) plotted at 
—6, —4, —2, 2, 4, 6, 8, 10 X a, where a is the noise level. The reference 
beams can be seen, due to the east-west chopping and nodding of the sec- 
ondary mirror on the JCMT. The positions of the components detected by 
our SMA imaging, discussed later, are shown as crosses. 



these two submm components, 16.6 ± 3.8 mjy, the equivalent of 
18. 7± 4.3 mJy at 850 fim (applying a +13-per-cent correction for a 
55-K greybody with emissivity index, +1.5). This is consistent with 
the 23.8 ± 3.5 mJy measured by SCUBA at 850 at the ~1-cj 
level. 

The submm morphology appears consistent with that seen at 
lower S/N using the Plateau de Bure Interferometer (PdBI) of the 
Institut de Radio Astronomie Mi llimetrique (IRAM) at 240 GHz 
(Fig. [4]- Papadopoul os et al . 2000) and confirms the extended na- 
ture of the dust emission. The measured size, even taking both com- 
ponents together, is significan tly smaller than the [11 x 4]-arcsec 
deconvolved size estimated bv lStevens et aljj2003l) . and the single- 
dish and interferometric flux densities would allow up to 5 mJy 
of emission distributed on scales larger than those probed by the 
SMA. Overall, our observations show l ittle compelling evidence for 
the galaxy-wide star burst reported bv lPapadopoulos et al J d2000f) or 
IStevensetal] ( l2003l) . 

3.2 CO morphology 

Papadopoul os et al] (|2000) found the CO J = 4 — 3 line emis- 
sion from 4C 60.07 to be bright and broad, even by the sta ndards 
now known for SMGs dNeri et alj|200ll ; iGreve et alj|2005l) . They 
describe two distinct CO(4-3) components in their relatively low- 
spatial-resolution PdBI data cube (Fig. |4j. Together, these cover 
>l,000kms -1 . 

One CO(4-3) component is spans only ~150kms -1 FWHM. 
Its position is consistent with that of the putative AGN core - taken 
to be the flat-spectrum radio source between the steeper-spectrum 
hot spots seen most clearly at 1.4 GHz - and it is -~220kms _1 
blueward of the systemic (He II) velocity of the AGN's host galaxy. 



Contours: SMA 890um 
10 Greyscale: IRAC 3.6+4.5um 





10 5 -5 -10 

ARC SEC 



Figure 3. Greyscale representation of the Spitzer averaged 3.6- and 4.5-/im 
images of 4C 60.07. Superimposed on the IRAC data is the high-resolution 
SMA 890-/xm image of 4C 60.07, with contours plotted at -2.5, 2.5, 3.5, 
4.5, 5.5 and 6.5 X 1.1 mjy beam -1 , the r.m.s. noise level. The size and 
shape of the SMA's synthesised beam is shown in the bottom left corner of 
this plot. The FWHM of the SMA's primary beam extends slightly beyond 
the extent of the image, to ~30 arcsec FWHM. The position of the radio 
core - as seen at 4.7 and 8.2 GHz - is shown by a square; those of the steep- 
spectrum hot spots are labelled with diamonds; the submm components, 'A' 
and 'B', are labelled. The IR, radio and submm images are aligned to better 
than 0.5 arcsec, conservatively. 



IGreve et "all J2004I) detected gas glowing in CO(1-0), coincident 
in velocity with the narrow CO(4-3) emission, and consistent spa- 
tially. On the basis of dynamical mass limits. IPapadopoulo s et al.l 
argued that this narrow component has yet to form the bulk of its 
eventual stellar mass. 

The second CO(4— 3) component lies -~700kms~ redward 
and ~7 arcsec south-west of the narrow component and spans 
^550kms~ FWHM, with emission li kely extending to lower fre - 
quencies, b e yond the observed band ^Papadopoulos et ai] |200rj|) . 
IGreve et alj ( 120041) report CO(1-0) in the velocity range of this 
broad CO(4-3) component, with a peak flux density of 0.27 ± 
0.05 mJy, but offset spatially by ~4 arcsec to the north-east, lying 
just west of the radio core. 

However, our re-analysis of the data used by IGreve etHll to 
explore the broad CO component, together with a much larger 
set of data taken with shorter baselines, have revealed no CO(l- 
0) emission above 0.05 mJy beam -1 (Fig. [5}. The total CO(1-0) 
flux, integrated over the velocity range to which these data are sen- 
sitive (+455 ± 312 km s -1 ), is 3<r < 0.07Jykms~ 1 . Adoptin g 
Xco = O.8(Kkms~ 1 pc" 1 )" 1 M toownes & Solomon|[l998h , 
this translates into Mh 2 < 3.6 X 10 10 M Q , less than half of 
the mas s determined on the basis of the CO J = 4 — 3 emis- 
sion by Papadopoulos et al. who assumed a 4—3/1-0 line ratio of 
r43 = 0.45, and less than the mass determined for the narrow 
CO(l- O) component seen blueward of our passband bv lGreve et al.l 
(2004). This suggests variations in the bulk gas properties across 
the 4C 60.07 system. 
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Figure 4. Greyscale representation of the uniform-weighted, velocity- 
integrated CO(4-3) emission from 4C 60.07, as detected using IRAM PdB I 
with a [8.9 X 5.5]-arcsec synthesised beam by Papadopoulos et al. (2000), 
superimposed with contours of their 240-GHz (1.25-mm) continuum im- 
age, plotted at -2.5, 2.5, 3.5, 4.5, 5.5 and 6.5 X 0.58mJybeam _1 , the 
r.m.s. noise level, with the shape and size of the somewhat elongated 240- 
GHz synthesised beam shown in the bottom left corner. The image is cen- 
tred on the position of the radio core - seen at 4.7 and 8.2 GHz - which is 
labelled with a square; those of the steep-spectrum hot spots are labelled 
with diamonds; submm components, 'A' and 'B' (see Fig. [5), are marked 
with crosses. 



3.3 Comparison of IR, submm and radio positions 

We first concern ourselves with the new data from the SMA. The 
most unexpected aspect of the SMA, Spitzer and VLA imaging 
(Fig. [3} is the offset between the primary sites of submm emis- 
sion and that of the radio core's svnchrotrorF^l The submm com- 
ponents lie 10-30 kpc (on the plane of the sky) from the radio core 
- a larger distance than can be understood via the positional un- 
certainties (<C4kpc; see Tableland ij2). The large offsets often 
seen between single-dish submm detections and their radio or IR 
counterparts dlvison etal, I l2007h are not expected here, due to the 
small SMA synthesised beam and accurate phase-referenced inter- 
ferometric positions. It is clear, therefore, that the radio galaxy is 
not associated directly with rest-frame far-IR emissio 

Fig. [6] shows how the morphology of the 4C 60.07 system 
changes as we move from rest-frame 0.2 (UV), via rest-frame 
0.5/zm(~ V),0.Sfj,m (~ R),0.9fj,m (~ /), 1.2 (~ J), 1.7 fim 



12 The positional offset is significant at the 4-<r level. The uncertainty in 
the submm positio n is dominated by the centroiding (see the appendix of 
Ivison et al. 2007) since systematic uncertainties are limited by the coin- 
cidence of the southern submm component with IR emission and by the 
positional tests described in i|2.1| 

13 Of the six submm-bright HzRGs with deep Spitzer imaging 
lArchibald et~aI1l200ll ; iReuland et aT1l2004l ; ISevmour et al]|2007l) , all dis- 
play an accurate alignment of the radio core and the 3.6-fj.m emission (to a 
tolerance of ~ 1 arcsec) and only 4C 60.07 shows a misalignment at 8 fim 
(Fig.|6). 
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Figure 5. Greyscale representation of the 24-GHz continuum emission from 
4C 60.07 superimposed with contours of the CO(1-0) emission in the ve- 
locity range of the b road CO(4-3) component detected with IRAM PdBI by 
IPapadopoulos et aU (+455 ±312 kms" 1 ) as imaged using the VLA with 
a [2.6 X 2.4]-arcsec synthesised beam (shown, bottom left). The CO(l— 
0) data shown are the average of the first two methods described in i|2.4l 
Contours are plotted at —2, 2 and 3 X 1 8.4 fjjy beam" 1 , the r.m.s. noise 
level. The emission coincident with component 'B' has an overall signifi- 
cance, using a Gaussian fit with size fixed to that of the synthesised beam, 
of -£,2 a. The image is centred on the position of the radio core, which is 
labelled with a square; those of the steep-spectrum hot spots are labelled 
with diamonds; submm components, 'A and 'B' (see Fig. [3), are marked 
with crosses. 



(~ H) and 3.3 pm (~ L) to rest-frame 5.0 (~ M), i.e. A b s 
= 0.9-24 /jm. The rest-frame UV emission is centred on the radio 
core. In rest-frame V, a red object ( 'K' ) becomes visible, 2.2 arcsec 
NNE of the core; the core emission covers ~2 arcsec and there is 
only low-level flux from the submm components. 'K' makes contri- 
butions in rest-frame R and / (Aobs = 3.6^1.5 fim) and the emission 
coincident with submm component 'B' strengthens. By the time we 
reach the rest-frame J and H filters at A D b s = 5.8-8.0 /im, 'B' has 
become the dominant component, with only very weak emission 
associated with the radio core and 'K'. Component 'A' appears to 
be devoid of significant mid-IR emission, with less than 4, 6, 26 and 
39/iIy (3<t) at A obs = 3.6, 4.5, 5.8 and 8.0 /^m. At A obs = 24 /jm 
we see emission along the orientation of the submm components, 
'A' and 'B', with little room for a significant contribution from the 
radio core (Fig. [6j. (The MIPS 24- filter transmits between 4.3 
and 5.4 /im - rest frame, to half pow er - and so misses the strongest 
PAH features, e.g. lDesai et alj2007l) . 

3.4 Clues from SEDs via IR diagnostic plots 

We have extracted I and 3.6-8.0-^m flux densities for the two 
components, the radio core and 'B', that dominate the rest-frame 
optical-IR emission f rom the 4C 60.07 sys tem, and adopted other 
measu rements from Chambers et al. ( 1996) and Ivan Breugel et al.l 
( 1998). The different noise levels at 3.6-8.0 in Fig. |6]make it 
difficult to judge the spectral energy distribution (SED) shapes, so 
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Figure 6. Greyscale representations of WHT / ( i|2.3h Keck K' I van Breugel et al. 1998), and Spitzer images of 4C 60.07 from 0.9 through to 24 /im, using 
a linear stretch from 0-10 <r, where a is the pixel noise (0-20 cr at 24 fim, where we have re-mapped the image into 0.86-arcsec pixels to match the IRAC 
images, 2.9 X smaller linearly than is necessary to fully sample the 24-fim point spread function). The WHT and Keck images have been smoothed with a 
0.3-arcsec Gaussian. The Keck image has been rotated slightly to match the USNO-based astrometry of the WHT /-band image. In each image, two crosses 
show the positions of the SMA submm peaks and the approximate rest-frame wavelength (in fj,m) is shown in parentheses. The images are centred on the 
position of the radio core, labelled with a square, while those of the steep-spectrum hot spots are labelled with diamonds. Each image is [25 X 25] arcsec. 



Table 3. Mid-IR flux densities. 



Isophotal wavelength 


Flux density (/^Jy) 


Qj,m) 


Core/A' 


'B' 


3.535 


20.6 ±1.4 


5.9±1.3 


4.502 


27.0 ±2.0 


12.1 ±2.0 


5.65 


3(7 < 26.0 


29.0 ±9.0 


7.74 


3cr < 39.0 


76.0 ± 13.0 


16.0 


175 ±33 


24.0 


930 ± 100 



we plot both of the dominant components in Fig.[7jand list the flux 
densities in Table [3] The rest-frame W-V emission from the core 
is resolved. At rest-frame R and / (A D b s = 3.6 — 4.5 /im) we see 
a blend of emission from the core and 'K'; since our photometry 
is appropriate for point sources, the true flux densities are probably 
~20 per cent higher than those shown. Component 'B' is barely re- 
solved (3.1 arcsec, FWHM), so its photometry ought to be accurate. 
The 24-/itn emission is a blend, with component 'B' contributing 
less than half. The totaJ3 flux density at 24 /jm is 0.97 ± 0.05 mjy. 

14 ISevmour et alj [2007) report different 8- and 24-fj.m flux densities for 
4C 60.07. Their 7-arcsec aperture, centred on the 3.6-fim emission, misses 
around half th e 8-/im flux (which comes mostly from 'B'). At 24 (im, 
ISevmour et alj used an aperture 26 arcsec in diameter (not 1 3 arcsec as 
stated in the paper; N. Seymour 2008, private communication) and thus 
include flux from the source to the south-east. Their 16-/xm aperture suffers 



AGN SEDs rise as a power law between 1 and 10 pirn 
(rest frame), while starbursts have flatter SEDs at rest-frame 1- 
3 /im which rise a t longer wavelengths (e.g. Ward et alj \m% 
iDonlev et ai]|2007h . At z = 3.8, IRAC probes rest-frame 0.8- 
1.7 ^m, so 5 , 8.o / jm/S'4.5/jm is a powerful AGN diagnostic, par- 
ticularly when p lotted against a probe of the >3 /xm slope, such 
as Sm um/Ssum dlvison et alj2004lLa"cv et alj2004l;lMagdis et al.1 
l2008h . 

If the 3.6-5.8-^m emission associated with the radio galaxy 
core is associated with starlight then we might expect that emis- 
sion to continue through to A b s = 7.7 /im - near the isophotal 
wavelength of the IRAC 8.0-/im filter, which corresponds to the 
l.6-[im stellar bump at z ~ 3.8. We do not see this, but using 
the IRAC photometry to constrain rest-frame //-band luminosities 
and assuming Lh/M ~ 2, the stellar mas s is poorly cons t rained : 
< 10 12 M , cf. (2.8 ± 1.2) x 10 11 M in lSevmour etai] ( 12007b . 
who found the //-band light to be 95-per-cent stellar. We note that 
a much larger L/M (and hence lower stellar mass) is plausible, via 
fine-tuning of age and initial mass function, and there is probably 
a contribution to the rest-frame near-IR light from the AGN. Our 
limit is therefore likely to be very conservative, but is consistent 
with a large mass of ~0.2-l-Gyr-old stars in the radio galaxy host 
accor ding to a simple Starburst99 model l lVazquez & LeithereJ 
2005). Already revealed as an AGN by its radio emission, the core 
does not have the optical-IR continuum characteristics of an AGN. 

In contrast, the SED of component 'B' climbs steeply to 



no such problem, although it was 12 arcsec in diameter rather than 6 arcsec, 
as stated, so we adopt their i6-fim flux here. 
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Figure 7. SEDs of the two main components of 4C 60.07 from the opti- 
cal to the submm, with component 'B' represented by open circles and 
the AGN core by filled circles. Beyond 8 fim we show total flux densi- 
ties, as diamonds. The dotted line is the S EP of the z = 1.44 ERO, 
HRIO telston et alj|l988t IStern et alj|2006l) . normalised in flux with re- 
spect to 4C 60.07 at rest-frame 185 fim (only the well-sampled optical/IR 
SEP is plotted); the dashed line is the SEP of a model 10 13 -Lq starburst 
lLagache et alj20Q3l) . 



Figure 8. Schematic of 4C 60.07. The horizontal axis is a spatial vector on 
the sky, starting at the radio core and running through 'A' and 'B'. The ver- 
tical axis shows relative velocity, whe re available. Black lines s hown the ex- 
tent of CO(4— 3) channel maps from Panadonoulo s et al .(2000) - the dashed 
line represents the emission that spans both components (note the velocity 
gap between it and the 150-kms -1 component). Components A' and 'B' 
are marked by circles at arbitrary veloc ities. The upper b ox represents the 
CO(1-0) map of the narrow component (Greve et al. 2004), while the lower 
box represents the 1.2-mm continuum, its velocity chosen to match that of 
the broad component. 



-^obs ~ 8.0 fim. It is fit well by a S„ oc v 33 power law: an 
extremely red galaxy, akin to HR 10 (= ERO J164502+4626.4) at 
z ~ 1.4 which has excellent co verage of its SED - see Fig. [7] 
jDev et alj[l999l : IStern et al.ll2006l) . The value of S 8 . ^ m /S 4 .5^ 
seen for component 'B' (6.3 ± 1.5, or [4.5] — [8.0] = 2.0) is un- 
usual : redder than any SMG in the SHADES sample ( Ivison et al. 
l2007h . and ranki ng with the r eddes t of the 8-/im-selected ob- 
jects explored by Magdis et al. I d2008l) . These diagnostics are typ- 
ically used for objects at z ~ 2, however, so we should explore 
whether the colour remains extreme at z ~ 4. Does it exhibit sim- 
ilar rest-frame optical/near-IR colours to HR10? For z — 1.4, 
S3.Gy.rn/ S2.iy.rn is close to S8.0Mm/S4.5 M m at z = 3.8. HR10 
has S3.6 M m/<S2.2 A Mn ~ 2.35 ± 0.24, so although HR 10 is redder 
than component 'B' at rest-frame ^1 /im, it is not nearly as red at 
rest-frame 1-10 fim which means that component 'B' can be cate- 
gorised as an obscured AGN with little room for ambiguity. 



4 A TOY MODEL OF THE 4C 60.07 SYSTEM 

It is tempting to invoke a complex, many-component system and 
to speculate that this form of interaction may be related to radio 
loudness and the formation of the most massive galaxies, with im- 
plications for the form and role of feedback. However, we must 
not forget the origins of our target: 4C 60.07 was selected initially 
as a USS emitter and its radio activity is thus inescapably young, 
^lOMyr dBlundell & Rawli ngs 1999). We are therefore seeing a 
system that is probably out of equilibrium. 

Can we marry together all the data for the 4C 60.07 system - in 
particular the radio and submm interferometry, both continuum and 
spectral line (CO), and the mid-IR imaging? These characteristics 
are presented schematically in Fig. [8] 

In our favoured model, 4C 60.07 comprises a far-IR-luminous 
starburst triggered by a galaxy-galaxy interaction. That more than 



one component is involved is inescapable, given the morphologies 
evident in the submm and mid-IR wavebands. Of the myriad com- 
ponents in the system, we identify the radio core as one with sig- 
nificant mass, since it appears to comprise a relatively large BH 
(capable of driving the radio lobes) and a host galaxy that is visi- 
ble in rest-frame UV, V, R and /. Since it lacks rest-frame far-IR 
emission, we suspect the radio galaxy host has exhausted (or driven 
away) its supply of molecular gas. Certainly, little H2 is evident via 
its most sensitive available tracers, dust and CO. We identify com- 
ponent 'B' as the other high-mass component, as evident from its 
very broad CO(4-3) line emission - a dusty, gas-rich galaxy, as ev- 
ident from the submm emission, undergoing a huge burst of star 
formation. Component 'B' is very heavily enshrouded in the dust 
that betrays its presence in the submm waveband, visible elsewhere 
in its SED only longward of rest-frame R, probably via a buried 
AGN with a power-law SED ( $3A[ . 

What of component 'A'? The key piece of evidence is its 
narrow CO line w i dth, s een in both J = 4 — 3 and 1 — 0. 
Papadopoulos et ai] feOOCh believed this emission was associated 
with the radio galaxy host but our high-resolution SMA imag- 
ing, together wit h the high-resolution CO(1-0) observations of 
iGreve et alj d2004l) . betray an association with component 'A' . We 
believe that 'A' is a bridge of cold material tidally stripped from the 
massive, gas-rich starburst/ AGN ('B') towards the host of the ra- 
dio core - an analogue of the scenario seen in some local ULIRGs, 
but on a more massive scale. This is in accord with its observed 
characteristics: the submm component has no counterpart at mid- 
IR wavelengths, is resolved on a scale of ~2arcsec (~15 kpc) and 
its gas - consistent spatially for CO(1-0) and CO(4—3) - has a very 
low velocity dispersion (~150kms _1 ). 
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5 CONCLUDING REMARKS 

New, high-resolution submm imaging of the HzRG, 4C 60.07, at 
z = 3.8, reveals two clumps of emission separated by 3.3 arcsec 
(~25 kpc). Although one clump is resolved by our ~2-arcsec 
beam, both are relatively compact (£15 kpc) and there is little evi- 
dence for a galaxy-wide starburst on larger scales from our obser- 
vations. 

Previously thought to comprise two gas-rich starbursts, one 
centred on the radio galaxy, we discover an extra level of complex- 
ity due to the misalignment of the submm emission and the radio 
core. Given the extreme youth of HzRGs, we interpret the 4C 60.07 
system as an early-stage merger, with its radio jets and a nearby 
gas-rich starburst/AGN both induced by the interaction. We inter- 
pret a second clump of dust, resolved by our submm imaging and 
coincident with very narrow CO emission, as a plume of cold, dusty 
gas in a tidal stream, seen in a non-equilibrium state. It is separated 
spatially from the radio galaxy, which may have exhausted its sup- 
ply of fuel. 

Recent work on the z = 4.7 quasar BR 12 02-0725 
dKlamer et all |2004|) . on LH 850.7 favison et alj l2002h and on 



Minkowski's Object dCroft et alj|200q) has suggested that some 
starbursts are induced by the interaction of a radio galaxy jet with a 
large gas reservoir. The story we are seeing unfold in the 4C 60.07 
system raises the possibility that some of these jets may have been 
triggered by the interaction that led to the starburst, rather than the 
starburst being triggered by the jet. 

When a BH is spatially offset from the gas and dust of its host 
system, as seen in 4C 60.07, this has implications for the calculation 
of its dynamical mass (via the extent of its CO, spatially and in 
velocity) and related claims ab out the co-assembly of B Hs and host 
spheroids at high redshift (e.g. lWalter et alj2003ll2004l) . regardless 
of whether the system in question is a disk or a merger. 

Our observations of 4C 60.07 bring t o mind scaled-up ver- 
sions of the NGC 4038/403 9 (Antennae' - I Wang et al]|2004l) and 
VV114 dFraveretalJl 19991) systems - early-stage mergers where 
the bulk of the star formation occurs in an obscured 'overlap re- 
gion' between two gas-rich nuclei. In NGC 4038/4039 this region 
displays a heavily enhanced star-formation rate yet contains lit- 
tle stellar mass (~10 per cent of the total - IWang et al.ll2004l) . In 
4C 60.07, the galaxy nuclei correspond to the radio core and 'B', 
and the overlap region is 'A'. There are also similarities between 
4C 60.07 and HE 0450-2958, one of the objects originally thought 
to have been caught during its transition from a ULIRG into an op- 
tically bright quasar |Huichings & Nefdl 19881 ; ISanders etalll988l ; 
ICanalizo & Stocktonkoo'll) but now believed to be a gas-rich spiral 
in collision with - and stimulating optical ac tivity in - an unob- 
scure d quasar hosted by a gas-poor elliptical (Papadopoul os et ail 
120081) . 

At z <; 3, where quasars are becoming rarer, a binary quasar 
with a small separatio n (< 1 arcmin) has yet t o be found (at 2 < z < 
3, tens are known -iMorflock et al .11 19991 ; iKochanek et all 1 19991 ; 
iHennawi et al.l2006h . The 4C 60.07 system contains two AGN and, 
although not strictly sp eaking a 'b i nary q uasar', it does support 
the argument made by IPiorgovskH dl99ll) that the high number 
of binary quasars on small angular separations - orders of mag- 
nitude more than that predicted by an extrapolation of the quasar 
correlation function pow er law - is due to the triggering of AGN 
during mergers (see also lSmail et~alll2003l ; [Alexander et al.ll2003l ; 
llono et al.ll2006h . 

We would not have been aware of the BH that presumably 
drives the immense radio luminosity of 4C 60.07 were it not for 



the interaction that triggered both the FRlI radio activity and the 
nearby starburst. Is there an unseen population of BHs at high red- 
shift - invisible unless actively fueled, yet sufficiently massive to 
drive a >10 27 -W radio source? Aside from its rare radio proper- 
ties, the SED of 4C 60.07's host galaxy is not particularly unusual, 
and such objects would be impossible to sift from the ~ 5 x 10 8 
simi larly bright IRAC sou rces in the sky. 

lReulandetalJd2007l) claim that HzRGs and SMGs differ only 
in the activity of their AGNs. Why, then, are complex submm mor- 
phologies the exception rather than the rule in typical SMGs? The 
answer lies in the youth of the 4C 60.07 system, as set by its selec- 
tion as a USS radio emitter. We see 4C 60.07 in the first throes of 
a violent interaction, probably within a few Myr of its latest burst 
of activity commencing, whereas analysis of the rest-frame opti- 
cal properties of SMGs s how that they are typically ^10x older 
than HzRGs - ~ 100 Myr dSmail et al.l2004l;ISwinbank et al.l2006h 
- and are therefore more likely to display more relaxed morpholo- 
gies. This provides a ready explanation for the exten ded nature 
of the submm emission seen towards the IStevens et alj sample of 
HzRGs , which is at odds with the compact emission seen f or most 
SMGs dYounger et ai]|2007l 120081 ; iBiggs & Ivisonl fcoOS). Since 
these HzRGs will all be youthful, we expect non-equilibrium mor- 
phologies to be commonplace, though not necessarily the complex 
dual-AGN signature seen in 4C 60.07. It is not the activity of their 
AGNs that mark HzRGs as different from SMGs, it is their youth. 
The possibility that HzRGs differ from SMGs has many implica- 
tions. It would be desirable to make similar observations of a larger 
sample of HzRGs at a resolution similar or better than that em- 
ployed here. 

Recent studies of rest-frame far-IR emission , either directly 
dTacconi et ail 120061 ; lYounger et al.l|200l |2008al) or via the far - 
IR/radio correlation dChapman et alj|2004 IBiggs & Ivisonll2008h . 
have highlighted the importance of hi gh spatial resoluti on. Single- 
dish observations - with SCU BA2 dHolland et al]|2006l), He rschel 
dGriffin et alj2007l) or SPICA dSwinvard & Nakagawal20o3) - will 
often provide an ambiguous or misleading picture. Our work re- 
inforces this view: we have spent a decade assuming that the radio 
core of 4C 60.07 dominates the rest-frame far-IR emission when the 
emission is, in fact, associated with a gas-rich companion. We must 
await sensitive interf erometers such as the Atacama Large Millime- 
tre Array (ALMA - IWoottenll2003h to probe the interactions that 
drive many aspects of galaxy evolution, and desig n future missions, 
such as the Far-Infrared Interferometer (FIRI -[Helmich & Ivisonl 
2008), with the ability to probe a large range of spatial scales. 

Facilities: SMA, Spitzer, VLA, ICMT, WHT 
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